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ABSTRACT 

The blazar 3C 279 is well studied and shows frequent large continuum flares 
from radio to 7-ray wavelengths. There have been a number of multi- wavelength 
observations of 3C 279, and hence there are many ultraviolet data for this 
object available in the UV archives. In this paper we present Lya emission line 
measurements for 3C 279 using all the archival lUE SWP spectra from 1988 to 
1996 and all archival HST/FOS G190H spectra from 1992 to 1996. 

Individual archival lUE spectra of 3C 279 show weak Lyo; emission at 
~1868 K = 0.536), which is easily seen in the co-added data. The Lya 
emission is observed in all the HST/FOS spectra. The strength of Lya is nearly 
constant (~5 xlO"^^ erg cm~^ s~^), while the 1750 A continuum varies by a 
factor of ~50, from ~ 0.6 to 31.6 xlO~^^ erg cm~^ s~^ A~^. 

The behavior of the Lya emission line flux and continuum flux is similar 
to that of the only other well observed blazar, 3C 273, which shows constant 
line flux while the continuum varies by a factor of ~3. This near-constancy 
of emission line flux in the two best-studied blazars, suggests that the highly 
variable beamed continuum is not a signiflcant source of photoionization for the 
gas. Some other source, such as thermal emission from an accretion disk, must 
be providing a significant fraction of the photoionizing flux in these objects. The 
large amplitude variability seen at 7-ray energies must be due to changes in the 
energetic electrons in the jet rather than changes in the external photon fleld. 

Subject Headings: galaxies: active - quasars: individual (3C 279) - 
ultraviolet :observation - variability 
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1. Introduction 

Blazars as a class of Active Galactic Nuclei (AGN) exhibit a smooth, highly variable, 
highly polarized IR-UV continuum. According to "Unified Schemes" for radio-loud AGN 
(Browne & Murphy 1987; Barthel 1989; Antonucci 1993) the observed continuum comprises 
a mixture of thermal radiation and non-thermal beamed radiation. The thermal, unbeamed 
component extending from the UV to the soft X-rays may originate from a hot accretion 
disk, as in radio-quiet AGN. The non-thermal component extending from the radio to the 
7-rays is thought to be relativistically beamed synchrotron radiation, possibly originating 
from a twin jet (Konigl 1981). This beamed component is highly variable and highly 
polarized. Blazars are thought to be those objects viewed at a small angle to the beaming 
axis and thus dominated by the non-thermal component (Urry & Padovani 1995). Some 
blazars exhibit a broad emission line spectrum similar to Seyfert galaxies and radio quiet 
quasars, particularly when the continuum is faint (the exceptions are classified as BL Lac 
objects; see however Scarpa & Falomo 1997). Although blazars are highly variable and show 
broad emission lines, their line variability studies have received relatively little attention to 
date. 

In recent years the study of broad-emission-line variability has been recognized as a 
powerful diagnostic of the broad line region in Seyfert 1 galaxies (Peterson 1993). In the 
best monitored AGN, high ionization lines like Lya and C IV are more rapidly variable 
(factors of ~10) than low ionization lines like Mg II or Ho; (e.g., Clavel et al. 1991; Peterson 
1993). Hence, the high ionization fines are extremely powerful probes of the smallest, most 
variable regions. Yet Lya emission line variability has been investigated in only a few 
high-luminosity radio-loud objects (Bregman et al. 1986; Perez et al. 1989; Gondhalekar 
1990; Ulrich, Courvoisier, & Wamsteker 1993; Webb et al. 1994; Scarpa, Falomo, & Pian 
1995, Wamsteker et al. 1997). Most of these studies indicate that emission line variability 
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timescales are on the order of months and the variabihty is confined to the fine cores. 
However, these objects are either not typical blazars (i.e., their continuum is not dominated 
by a beamed component) or have only a handful of observations. 

It is not clear whether the broad emission lines in blazars behave differently from the 
broad emission lines observed in Seyfert 1 galaxies. Further, what is the role of the blazar 
jet in photoionizing the emission lines; that is, what is the effect of the beamed blazar light 
on the line emitting region? What is the structure of the broad line region gas? Some 
investigations are starting to address these questions (e.g., Corbett et al. 1996). 

3C 279 is an excellent target to study emission line variability in blazars, since it is 
a highly variable, bright radio source. It is also one of the brightest known extragalactic 
7-ray sources (Hartman et al. 1992; von Montigny et al. 1995), is highly variable at X-ray 
and 7-ray energies (Makino et al. 1989; Kniffen et al. 1993; Maraschi et al. 1994), and is 
one of the clearest cases of relativistic beaming (Dondi & Ghisellini 1995). It shows highly 
variable optical polarization (9-45%; Angel & Stockman 1980; Impey & Tapia 1990; Mead 
et al. 1990) and has been one of the most intensively monitored objects of the blazar class. 
Finally, it is one of two emission-line blazars with extensive archival UV spectra. In this 
paper wc report the Lyct emission line history for 3C 279, which has not been as thoroughly 
investigated as the continuum, and its relation to the continuum variability. The analysis is 
based on archival ultraviolet data from both the International Ultraviolet Explorer satellite 
(lUE) and the Hubble Space Telescope {HST). 



2. Observations 
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2.1. HST/FOS Data 

Between 1992 and 1996, 3C 279 was observed three times with the HST Faint Object 
Spectrograph (FOS). The 1992 observation covers the wavelength range 1150 A - 4750 A, 
and has been discussed in great detail by Netzer et al. (1994) who note that the spectrum 
shows the Lya, C IV A1549, C III] A1909, Mg II A2800, and H/3 lines. The Lya, C III], 
and H/3 line profiles are fairly symmetric, while C IV and Mg II have substantial red wings. 
During this period 3C 279 was in a faint state, with a clearly visible Lya emission line (see 
Figure 1). The 1994 and 1996 observations cover the wavelength range 1150 A - 2312 A. 
Both these observations were obtained when 3C 279 showed a flare. Once again, in these 
FOS spectra the Lya is still easily visible (see Figure 1), but the signal-to-noise ratio (S/N) 
is not as good as in the 1992 observation. 

We retrieved the FOS data from the HST archive and recahbrated them using the 
best reference files, as recommended by the HST Data Handbook (Version 2). The details 
of each observation are given in Table 1 and the spectra are shown in Figure 1. The 
recalibration of the FOS spectra was necessary because the older spectra are cahbrated 
using different procedures and could not otherwise be directly compared to more recent 
spectra. The recalibrated FOS spectra can also be directly compared with the lUE spectra, 
since both the HST and lUE spectro-photometric calibrations use the model white dwarf 
reference scale (Colina and Bohlin 1994). In fact, based on an analysis of Seyfert galaxies 
observed simultaneously with FOS and lUE, absolute flux differences are ^ 5% (Koratkar 
et al. 1997). 

The 1994 FOS observation was obtained using two gratings (G130H and G190H). 
During the latter observation, due to reacquisition problems, the target was partially out 
of the aperture and the G190H spectrum was affected. We have therefore added a constant 
flux to the G190H observation to align the data with the G130H spectrum. The errors 
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given in Table 1 for this spectrum (y2ey0606) reflect our uncertainty (10%) in applying this 
procedure. 

2.2. lUE Data 

Between 1988 and 1996, 3C 279 was observed with lUE many times with both the 
short- wavelength and long- wavelength cameras. The long wavelength (LWP/LWR) spectra 
are affected by low detector sensitivity, and, increasingly with time, by scattered light, hence 
emission lines longward of the Lya line are difficult to detect. In this paper, we concentrate 
on the Lya line, and thus only the 15 short-wavelength (SWP) camera observations are 
discussed further. 

A careful inspection of the 15 line-by-line spectral images for flaws showed that in 
one case (SWP 40489), the Lya spectral region was contaminated by a cosmic ray. This 
spectrum was therefore rejected from further analysis. 

All the line-by-line spectral images were also analyzed to determine the background 
contribution. A comparison of the average Flux Numbers (FN) in the 1866-1872 A interval 
centered on the object spectrum and on the background (10 lines away on either side of the 
object spectrum) showed that in SWP 46662 and SWP 53261 the background contribution 
exceeded 90% of the signal-|-background emission. Therefore the Lya emission hne flux 
measurements for these spectra in Table 1 are 5cr upper limits. These noisy spectra have 
not been used in the co-addition described below. 

All the spectra were extracted from the line-by-hne image flies using the NEWSIPS 
technique being implemented in the flnal lUE archive. A summary of the NEWSIPS 
extraction technique can be found in Nichols & Linsky (1996). The details of each 
observation are given in Table 1. All the 1992-1993 spectra are under exposed, but some 
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signal is present. lUE spectral regions with known artifacts, at 1278-1279 A, 1450-1550 A, 
and 1663 A (Crenshaw, Norman, & Bruegman 1990), were excluded from the analysis. 

The improved S/N obtained through NEWSIPS allowed sufficient contrast in the SWP 
spectra between the continuum and Lya, even in the lowest states, to perform for the first 
time a systematic study of the emission line in this source. The Lya emission feature that 
is clearly present in the FOS data is also seen in all but four individual SWP spectra (see 
Table 1). Recently, the sensitivity degradation correction for SWP spectra obtained after 
1993 has been determined. This correction is time dependent, and can be up to 5% in 
1996 (Imhoff 1997). This correction has not yet been implemented in the final lUE archive, 
however, so we have conservatively added a 5% error in quadrature to the uncertainties of 
the continuum and hne fluxes for spectra obtained after 1993. 

To improve the S/N ratio of the Lya emission line, we binned the spectra in four 
groups depending on the continuum flux at 1750 A. The SWP spectra in each group were 
then co-added, as follows: the 'high spectrum' was generated by co-adding all spectra with 
fcont > 2 X 10~^^ erg cm~^ s~^ A~^, the 'medium spectrum' is a co-addition of all spectra 
with 8 X 10~^^ < fcont < 2 X 10~^^ erg cm~^ s~^ A~^, the 'low spectrum' is a co-addition 
of all spectra with 3 x 10^^^ < fcont < 8 x 10~^^ erg cm~^ s^^ A^^, and the 'very low 
spectrum' includes all spectra with 1 x 10~^^ < fcont < 3 x 10~^^ erg cm~^ s~^ A~^. In 
Figure 2 we show the four co-added spectra. The Lya emission line is easily seen at ~ 1868 
A. The redshift measured from the lUE observations is in good agreement with the redshift 
determined from the high S/N ratio FOS data (Netzer et al. 1994) given the wavelength 
uncertainty and the low S/N ratio of the lUE spectra. 



-8- 



2.3. Flux Measurements 

The Galactic hydrogen column density in the direction of 3C 279 is 2.22xl0^°cm~^ 
(Elvis, Lockman, & Wilkes 1989) which implies Ay = 0.13 for a total-to-selective extinction 
ratio R = 3.1. All the FOS and lUE spectra were dereddened assuming the Cardelli, 
Clayton, & Mathis (1989) extinction curve. 

A simple power-law model was fitted to the dereddened spectra in the wavelength 
region 1230 - 1950 A using the IRAF/STSDAS tool nfitld. Wavelength regions affected by 
the lUE artifacts were excluded from the fits to the lUE spectra, as was the region around 
Lya. We defined 1750 A as the fiducial wavelength for continuum fiux measurements 
because this region of the SWP spectrum is relatively free of lUE artifacts and line emission. 

After subtracting the best-fit power-law continuum, the Lya line was measured as 
follows. The 1992 observation, which has the best S/N and spectral resolution, was 
accurately modeled with a multicomponent Gaussian fit using a minimalist approach (least 
number of free parameters). IRAF/STSDAS tool ngaussfit was used for this procedure. 
This model of the Lyo; line was then used to fit the remaining two FOS spectra. Since 
these spectra have lower S/N, the observed wavelength of the Lyo; line was fixed at 1868A 
(within the wavelength calibration error of the FOS) and the widths and amplitudes were 
allowed to vary. We note here that all the FOS spectra show an intrinsic absorption feature 
in the Lyo; emission line. These absorption features are discussed in detail by Chaffee et al. 
(1997). 

Due to the low S/N ratio of the lUE data, the observed wavelength of the Lya line 
was fixed at 1868 A (within the wavelength cahbration error of the lUE), and the line was 
modeled with a single Gaussian whose width and height were allowed to vary freely. A 
single Gaussian fit was sufficient for the lUE because of the low spectral resolution. In all 
but two spectra, the Lycu emission fine is easily detected; in the exceptions (SWP35443, 
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SWP36420) the Lya region is very noisy. 

The continuum fluxes at 1750 A and the Lya line fluxes, along with their associated 
uncertainties, are reported in Table 1 for the individual spectra and in Table 2 for the 
co-added spectra. 

3. Results and Discussion 

Although the data arc sparse and unevenly sampled, the continuum flux density 
and the line flux over the 8-year observation period show very different characteristics 
(Figures 3 and 4). The continuum flux has varied by a factor of ~50, while the Lya line 
flux has remained nearly constant. We used two independent methods to test whether 
the emission line flux is constant. In the flrst method, a reduced chi-square test on the 
emission line flux gives a = 1.3 for 12 degrees of freedom indicating that the emission 
line flux data are consistent with being constant. In the second method, in our line fltting 
procedure (discussed in section 2) we kept all parameters of the Gaussian model constant 
and determined the reduced chi-square of the flt for each spectrum. The constant model 
used the parameters derived from the 1992 FOS spectrum, which has the best S/N and 
spectral resolution. The reduced chi-square of the constant model and the free Gaussian 
model were checked using an F-test. This test showed that the constant model and the 
free Gaussian model were comparable at the 95% level (F=l.l for a change of 2 degrees of 
freedom). We conclude that in 3C 279, although the observed continuum flux density has 
varied, the Lya flux has remained essentially constant. 

This trend between the continuum and line is similar to the case of 3C 273, where the 
Lya varied little or not at all during a factor of three change in the UV continuum (Ulrich 
et al. 1993). Such a trend between the emission line fluxes and the continuum was flrst 
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noted by Sandage, Westphal and Strittmatter (1966). The continuum and line variation in 
3C 279 is in sharp contrast to that seen in Seyfert galaxies, where the continuum and line 
vary with similar amplitudes and in a correlated way (Koratkar & Gaskell 1991, Clavel et al. 
1991, Peterson 1993, Reichert et al. 1994, Crenshaw et al. 1996 and Wanders et al. 1997). 
This lack of Lyo; variation is not just an effect of temporal sampling of the variability. In 
order to determine the effect of temporal sampling on the line variability, we compared our 
results to those of the much better-understood case of the Seyfert galaxy NGC 5548 which 
shows a clear correlation (r = 0.6, p ~ 99%) between Lya and the UV continuum (Clavel 
et al. 1991; Korista et al. 1995), with a lag of ~10 days. In 1989, NGC 5548 was observed 
every four days for eight months (Clavel et al. 1991). We selected 16 points at random from 
the 1989 NGC 5548 continuum and Lya light curves. For these 16 points we determined 
the amount of variation in both the continuum and lines flux. The process was repeated 
a 100 times. There was not a single 16 point sampling where the line flux was constant 
while the continuum flux varied as seen in 3C 279. This difference in the Lya variability 
of blazars can be explained if the thermal continuum that primarily photoionizes the line 
emitting gas is hidden by the variability of the beamed continuum in radio-loud AGN. This 
effect would be strongest in blazars, where the beamed component contribution dominates 
the observed continuum. 

The Lyct equivalent width (W^Lya) of 3C 279 ranges between lA and 45A (Figure 
5). These values of WLya are smaller than those observed in normal quasars (Kinney, 
Rivolo & Koratkar 1990; Cristiani & Vio 1990; Cheng & Fang 1987), where the thermal 
component, possibly the accretion disk, is the only source of ionizing radiation and observed 
continuum and gives rise to the observed Wiya ~100A. In blazars there are two sources 
of ionizing radiation: (1) the beamed non-thermal radiation from the relativistic jet, 
and (2) the thermal component probably from an accretion disk. Thus, the observed 
continuum has contributions both from the disk and jet. If we assume that the beam 
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ls not contributing to ionizing the line emitting gas, and only the thermal component 
is responsible for line emission, then the contribution of the beamed component to the 
observed continuum determines the strength of the equivalent width. The expected Wiya 
in blazars should be smaller than observed in normal quasars and should be anti-correlated 
with the contribution of the beamed component to the observed continuum (see Corbett et 
al. 1996). In 3C 279 we see this characteristic equivalent width trend (see Figure 5). At 
low observed continuum levels, the WLya becomes comparable (~100A) to other quasars 
when the beamed continuum contribution is ~2x the accretion disk contribution, and at 
high observed continuum levels the beamed continuum contribution is ~30x the accretion 
disk contribution. This model suggests that at low continuum levels we should see a change 
in UV continuum slope because of the increased contribution of the big blue bump due to 
the accretion disk. At high continuum levels we should just see the synchrotron power law 
slope. Such a change in the UV slope of 3C 279 has been noted by Pian et al. (1997). 

One of the models for the production of 7-ray emission in blazars is Compton 
upscattering of ambient Lya line photons (Sikora, Begelman, & Rccs 1994). Variability 
of 7-rays could therefore be caused by changes in the Lya intensity or in the energetic 
electrons of the jet. The 7-rays would vary linearly with Lya flux and/or quadratically 
with respect to synchrotron emission depending on the direction of the jet electrons. If 
some of the Lya flux increases in response to the beamed continuum (Ghisellini & Madau 
1996), then the 7-ray flux could increase even more than the square of the observed UV 
synchrotron flux. A strong correlation between UV and 7-ray flux has been observed in 
3C 279 (Maraschi et al. 1994, Wehrle et al. 1997), indicating that the beamed continuum 
could be ionizing some of the line emitting gas. Our present analysis suggests that 7-ray 
variability must be caused by variations in the jet rather than the external Lya photons. 
The present data do not allow us to determine the contribution of jet synchrotron emission 
to photoionizing some of the Lya emission clouds. If the beamed continuum does ionize 
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some of the line emitting gas, we should see rapid variations in the emission line profile. 
The nature of the variations would depend on the kinematic structure of the line emitting 
gas. Our lUE data do not have high enough spectral resolution or S/N to allow us to 
decompose the line core and wings but the Lya line profile appears symmetric, as it does 
in the FOS 1992 spectrum (Netzer et al. 1994). The 1994 and 1996 FOS spectra show 
intrinsic absorption and do not have sufficient S/N to determine accurately the line profile 
variations. Simultaneous monitoring of 3C 279 in both the UV and 7-ray would allow us 
to study continuum and emission line profile variability, and then unambiguously select the 
correct inverse Compton model for 7-ray production. 

4. Conclusions 

Our analysis of the archival lUE and HST/FOS observations of 3C 279 over a period 
of eight years shows that although the continuum has varied by large amplitude, the Lya 
fiux has remained essentially constant. This suggests that the beamed continuum is not the 
dominant source of photoionization of the broad line clouds, which are instead likely to be 
powered by a thermal component underlying the beamed continuum. 

The strong 7-ray emission from 3C 279 may be due to Compton scattering of the 
Lya line photons. In this case, the large amplitude of the 7-ray variability must be due 
to changes in the energy distribution of electrons in the jet, rather than changes in the 
external photon field. This implies that 7-rays will change as the square of the synchrotron 
flux, or perhaps even more if the beamed continuum ionizes line-emitting gas lying close to 
the jet axis (Ghisellini & Madau 1996). High quality line profile observations are necessary 
to test this model further. 
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Table 1: Fits to the Individual Spectra 



Image/ 


Date 


Exp. time 


F\ 

^1750 




Spectrum 


dd/mm/yy 


(sec) 


(xlO^-*^^ erg cm~^ s^^A^^) 


(xlO^^'^ erg cm^^ s^-*^) 


SWP33864 


05/07/88 


6000 


29.43±0.36 


12.17±4.96 


SWP33865 


05/07/88 


10200 


31.55±0.29 


6.63±1.95 


SWP35443 


28/01/89 

/ / 


10740 


13.83±0.21 


11.20 ^ 


SWP36420 


08/06/89 


13200 


17.09±0.17 


16.35 ^ 


SWP40489 


29/12/90 


14400 


18.34±0.14 


. . . b 


SWP42132 


27/07/91 


14400 


11.69±0.11 


8.09±1.60 


SWP44806 


29/05/92 

/ / 


10800 


12.24±0.26 


11.19±3.50 


SWP46649 


02/01/93 

/ / 


13200 


1.42±0.12 


4.52±1.29 


SWP46653 


03/01/93 


36000 


1.47±0.11 


3.88±0.98 


SWP46657 


04/01/93 


39600 


1.27±0.10 


5.33±0.70 


SWP46662 


05/01/93 


39600 


1.03±0.10 


3.72^ 


SWP49681 


24/12/93 


12000 


3.45±0.21 


7.38±1.33 


SWP49686 


25/12/93 


19800 


4.19±0.24 


4.44±1.15 


SWP53261 


03/01/95 


16800 


0.60±0.10 


5.38 ^ 


SWP56635 


25/01/96 


19500 


4.59±0.31 


5.27±1.55 


y0pe0b02/3 


08/04/92 


2250 


1.89±0.03 


5.22±0.56 


y2ey0606 


05/06/94 


1980 


6.12±0.08 


5.86±1.41 


y32j0605 


08/01/96 


2250 


5.71±0.03 


4.04±1.19 



^ The amplitude of Lya is less than the rms fluctuations about the best-flt power-law model 
so values given are Scr upper limits. High background seen in some of the spectra. 
^ Cosmic ray hit exactly on the Lya line, flux not measured. 
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Table 2: Fits to the Co- Added Spectra 



Image 


Co-added spectra 


(xlO~^^ erg cm~^ s""^ A"^) 


(xlO~^^ erg cm~^ s~^) 


'HIGH' 


SWP33864, SWP33865 


30.66 ± 0.47 


10.86 ± 6.45 


'MEDIUM' 


SWP35443, SWP36420 
SWP42132, SWP44806 


13.71 ± 0.36 


5.00 ± 2.45 


'LOW 


SWP49681, SWP49686 
SWP56635 


4.12 ± 0.14 


5.26 ± 2.12 


'VERY LOW 


SWP46649, SWP46653 
SWP46657 


1.40 ± 0.11 


4.20 ± 2.11 
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Fig. 1.— The three HST/FOS spectra of 3C 279. Note that even in the low state, the Lya 
hne is clearly visible at 1868 A. (a) 1992; (b) 1994; (c) 1996. 

Fig. 2. — The co-added lUE spectra of 3C 279 in four intensity states, as defined in § 2.2 
(see also Table 2). In all continuum flux ranges the hya line is visible at 1868 A. (a) 'High' 
spectrum; (b) 'medium' spectrum; (c) 'low' spectrum; (d) 'very low' spectrum. 

Fig. 3. — (a) The continuum flux density (filled points) and Lya line flux (open points) 
of 3C 279 as a function of time. Over the eight years of observations the source faded by 
a factor of ~ 50, followed by a slight increase at the most recent epoch. The lUE data 
are represented by circles and the HST/FOS data by squares. Error bars arc smaller than 
the symbol if not indicated in the figure, (b) Same as (a) but for co-added lUE spectra 
(pentagons) and HST/FOS data, (squares). 

Fig. 4. — (a) The Lya line flux of 3C 279 as a function of the continuum flux density (circles, 
lUE; squares, HST/FOS). (b) The Lya line flux of 3C 279 as a function of the continuum 
flux density for the co-added spectra and the HST/FOS (pentagons, lUE coadded; squares, 
HST/FOS) Over a period of 8 years the continuum flux varies by a factor ~ 50 while the 
line flux remained nearly constant. 

Fig. 5. — The Lya equivalent width as a function of the continuum flux density (circles, 
lUE individual; pentagons, lUE coadded;squares, HST/FOS). Such a decline in equivalent 
width is possible if the line is primarily ionized by the thermal component, but the flux in 
the observed continuum has varying contributions from the beamed continuum. The solid 
line represents the variation in equivalent width for a constant line intensity and a varying 
beamed continuum. 
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